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We obtain a simple formula for the semiclassical series for potentials V(x) = x N (N even) 
and derive almost explicit formulae for the WKB approximation of the energy eigenvalues 
of such potentials. 

In the present paper we study the structure of the WKB series for the polynomial 
potential V(x) = x N (N even). In particular, we obtain relatively simple recurrence 
formula of the coefficients a' k of the semiclassical approximation and of the WKB 
terms for the energy eigenvalues. We follow our aproach as expounded in ref.0'. 

We consider the two-turning point eigenvalue problem for the one-dimensional 
Schrodinger equation 

d 2 

[-?i 2 ^2 + V(x)]iP(x) = Eip{x). (1) 
We can always write the wavefunction as 

ij(x) = exp(ia(x)j (2) 



h 

where the phase cr(x) is a complex function that satisfies the differential equation 

a' 2 (x) + ha"(x) = (V(x) - E) d = f Q(x). (3) 
The WKB expansion for the phase is 

oo 

a(x) = Y / h k o- k (x). (4) 

k=0 

Substituting (HI) into (||) and comparing like powers of % gives the recursion relation 



'2 



Q(x), Vk^n-k + °n-l = 0- (5) 

k=0 



We assume that the potential V(x) is single valued analytic function which is 
real on the real axis, that V(±oo) = oo, that V{x) has a unique minimum on the 
real axis, that V{x) rises monotonically on both sides of the minimum and for k > 1 
the integrals 

/a /"+00 
\o-' k (x)\dx, / \a' k {x)\dx (6) 
-oo J b 
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are convergent, when a < x%,b > X2, where x^< £2 are the turning points. Then 
the following asymptotic expansion takes placeU: 



1 f 00 

— ® ^2 fi k a' k (x)dx = nnh, (7) 

2i Jy , 



2i ^ fc=0 

where n > is an integer number and 7 is a contour surrounding the turning points 
on the real axis (for the first time this formula appears, apparently, in |3|)). This 
relation is an equation with respect to E and using it one can find the asymptotic of 
the eigenvalues E n (K). In some cases the series (|7|) can be summed exactly (see e.g. 

!),§,§)J>!))- 

The zero-order term of the expansion (||) is given by 

^Jda = J dx^E-V(x), (8) 

the first odd term is 

and to find the higher order terms we need to compute the functions a' k using the 
recursion relation (||). 

We now consider the potentials V(x) = x , where iV is even positive integer. Let 
us show that for such potentials the coefficients a' k (k > 1) of the WKB expansion 
have the form 



a' 



(_z\3k-l„-k+N fc-1 
1 S . ; ^ , : X -h , :,/: / " ' : r'\ (10) 



(E-x 



1 



3=0 



where we choose VE — x N = iV x N — E, and the coefficients A k -j-ij of the mono- 
mials 3 ' -1 :e j are computed according to the recurrence formula 

^ s 2-1 

A s,l = o E E A i,i A s-i,l-l-j+ (11) 
i=0 j=0 

1(2 + N) + (2 + 3N)s - N (N-l)l + N-s. 

T As,Z-l H ^ A-s-1,1- 

We prove formula (11) by induction. We have 

a' = iVE^, *'x = A{E _ xNy (12) 

it means that for k = 1 formula (11) takes place, and also Aq^q = N/4. Let us 
suppose that it holds for k < m. Then for k = m in the case < u < m we have 

( -\3m-2 ~m+2N u—l m—u-l 



E A Tpu—j — ^rrjN \ ^ a rpm—u—r — 1 

s± u —j — ljr^ X / ^ 1 m-ii-r-l,r- c ' 

(13) 
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Note that for a fixed i,j the coefficient Aij is a coefficient in the expresion (|iTj|) for 
a' t if and only if i + j = t — 1. Then from (|l3| ) we get that the contribution into 



{e - x N y 

(s + I = m — 1) arises from 



' (14) 



_ 1 \3m-2 x -m+2iV 



E — : ^^A J -^^_ itI _ i _ 1 ^- < a; ^- 1 ) JV = (15) 

i+i=tt-l (£-a^)~ 



(_i)3m-2 x -m+A r 

3m-2 /, -AiJ-As— 1,/— J — 1-^ ^ ) 

(E-x N )s~ i+j =u-i 



here s — i + Z — j — 1 = m — u — 1 and, therefore, j— i is a coefficient in cx^„ u . 

Next, 

( _-\3m-4 -m+N m-2 
= -— a^T- EKl-^ + ^V + mAm-j-ljV"- J '^ jN + (16) 

(-1 - 3iV + m + ^miV - jiV)A m _ J -_ 2 j.E m - J '- 2 x w+1 ) JV ]. 

Hence, from (||), (|l~5D , (p^) we conclude, that A s j is computed according (11) and 
a' m has the form (|i0|). 

Consider the formal series 

A(zi,z 2 )= E ^,i44 (17) 

where Aij are numbers defined by the recurrence formula (11). Using the properties 
of generating functions (see, e.g.EP) we obtain from the recurrence relation (11) the 
following differential equation for the generating function A(z\,z 2 ) 



1 ,2 2 + N / am 2 + 3iV 

-z 2 A + -^—z 2 {z 2 A) Z2 + — - — z lZ2 A Zl 



N N - 1 1 N N 

-z 2 A + -^— Zl z 2 A' Z2 - -z x {z x A)' Zi + -z x A + — . (18) 

We cannot solve this equation in the general case, but we can use it in order to 
find the coefficients of the form A^o- Indeed, if we put z 2 = 0, then equation (18) 
has the form 

z 2 A'(z) - (z{N - 1) - 2)A{z) - N/2 = 0, (19) 
where we set z = z\. The general solution of this equation is 

A(z) = CeH 1 -^ lo sW + 2~ l ~ N e^Nz"- 1 r(N, -\ (20) 
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where r(n,z) = f z °° t n ~ 1 e~ t dt. However to satisfy the initial condition A(0) = iV/4 
we have to put C = 0. Therefore the coefficients A^o are the ones of Taylor expansion 
of the function 



"iVgf AT JV-1 2 



N z N ~ l r(N, -). (21) 



Integrating we get for integers N > 

2 J 

-) = L 2 ^-fc+2(A; - 1 



2, A AT! 



2-1-" ei N z^ r(N, ±) = £ oJV _ fc+ ^_ lM ^- fe . (22) 



Therefore 

N\ 

As >° = 2 s+2(N- s -iy: (23) 

and this is in agreement with ^4o,o = N/4. Unfortunately, we are not able to find 
the solution of the general equation (18) and, hence, cannot get similar expressions 
for other coefficients A s ^. 

However using relations @, (11) we can obtain the following formula for the 
coefficients of the WKB expansion 

i2 3fc + 1 ^F r{^ + l) 2 t r \3-2k 

/^ dx = W^WN E N W^TW) { Ajko S(— + 

I -2k . 2 ^ t J-. 1 - 2k 2k T t 1 ,3 - 2k 1 - 2k 

— - 8 )+x;^-i-i l in( a +— ) n (—J— +^v^" s) ' 

i=l s=l s=l 



(24) 



/From (M) we conclude that the integrals (pi) are convergent for all k > 1 and 



where /c > 1 and A<ik-i-\,i are computed according to (11) and 

2i£5+^0Fr(l + i) 



according to the above mentioned statement fromEP formulae (^) and (24) can be 
used to find the asymptotic approximation of the eigenvalues E n (K). We should 
emphasize that quite generally the odd-order terms (except first order for a[) yield 
integrals that vanish exactly, because, as it follows from the results of FromanlilP, 
a 2k+i are total derivatives for all k > 0. 

To summarize, in this paper we derive relatively simple explicit formula for the 
WKB terms for the energy eigenvalues of the homogeneous power law potentials 
V{x) = x , where N is even. In computing by means of this formula we manipulate 
only with numbers and do not need to collect similar terms of a polynomial, which 
we must do otherwise when we use just the recursion formula (|5|). As it is known 
the operation of collection of similar terms is very laborious even for Computer 
Algebra systems. So, application of the obtained formulae considerably simplifies 
calculations, especially if we need to compute high order terms. 
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